Diffusion induced mode switching in avalanching systems 
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Most avalanching systems in nature should involve diffusive processes as well which can change the behavior 
of such systems and should be taken into account. We examine the effects of diffusion on the model of a 
dissipative bi-directional burning model. It is shown that the system behavior progressively changes with the 
increase of diffusion. Avalanches of small sizes become suppressed and large quasiperiodic bursts develop. In 
contrast with sandpiles these bursts are not edge triggered and the transition is gradual. 
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In the last fifteen years it has been suggested that many 
complex physical systems behave like avalanching systems 
and follow the Self Organized Criticality (SOC) paradigm 
JlJ. This concept proposes a new statistical mechanics for 
open many-body systems with non linear dynamics, imply- 
ing that such systems spend most of the time near the crit- 
ical state. The SOC state if often characterized as the dy- 
namical state of an avalanching system in the limit of weak 
driving where a number of statistical parameters are char- 
acterized by power law distributions. SOC has been sug- 
gested to occur in a large number of very different physical 
01301113131111 and non-physical EH systems 
and phenomena. Among other systems where avalanching 
behavior is established or suspected, are a number of space 
plasma systems and phenomena like solar flares and current 
sheet. In these systems avalanches are related to the local- 
ized magnetic reconnection. There is growing understanding 
that avalanches may be a common feature or a rather wide 
class of systems. Most of the studies of such systems concen- 
trate on the avalanche statistics and possible transition to SOC, 
without invoking additional physical mechanisms of trans- 
port. However, there is a strong basis to believe that system 
dynamics can be changed in the presence of physical mecha- 
nisms like diffusion. Most of the models that were suggested 
until now hardly dealt with the influence of diffusion on the 
dynamics of the avalanche. Yet observations show that vir- 
tually all avalanching systems in nature involve processes of 
diffusion, weak or strong. Such observations lead us to believe 
that diffusion is fundamental process in any study of dynamics 
of these systems. The first attempt to add a diffusion flux in 
avalanching system was suggested in Ref.[3 The study was 
based on the one-directional discrete sandpile model. It was 
shown that the system breaks into two spatially distinct re- 
gions: diffusion dominated and avalanche dominated. For suf- 
ficiently strong diffusion the system switched into the regime 
of quasiperiodic large bursts, which were claimed to be edge- 
triggerred. Preferential direction governs the dynamics of the 
sandpile and makes the conclusions hardly applicable to other, 
isotropic, systems. In the present paper we study the influence 
of diffusion on the avalanches in a continuous bi-directional 
model. We show that diffusion does not break the system spa- 
tially into avalanche dominated and diffusion dominated parts. 
However, competition between diffusion and driving can force 



the system to switch from the regime of usual avalanching 
dynamics, where avalanches of various sizes are present, to 
the regime of large bursts with almost no small avalanches. 
The transition is not discontinuous but occurs gradually when 
avalanches of smaller sizes are suppressed as the diffusion be- 
comes stronger. 

We consider the one-dimensional bi-directional burning 
model 113LI14II . which was proposed as a more accurate (al- 
beit phenomenological) description of locally reconnecting 
plasma systems. The model is formulated in terms of energy 
transfer but actually was aimed to represent current carrying 
particles in the current sheet. The main variable is the temper- 
ature field. When the local temperature T(x,t) exceeds the 
upper critical value, T > T c , burning starts locally, during 
which energy releases at the rate Q = kT. Part of this en- 
ergy a is isotropically transferred to the neighbors while the 
remaining part (1 — a) is radiated out of the system. This 
radiated energy flux represents direct energy losses from the 
system and it is what can be measured by a remote observer. 
Because of the energy losses (flux to neighbors and radiation) 
the local temperature decreases until it gets below the lower 
critical value, 7], and the burning ceases. Thus, the smallest 
single site avalanche lasts for r m i„ w ln(T c /T;)/fc. Burning 
can start again only when the temperature increases beyond 
the upper threshold. This hysteresis behavior is the necessary 
feature of avalanching systems, here it represents the fact that 
the threshold for instability (temperature for the burning igni- 
tion or current for the reconnection triggering) is higher than 
what is needed for energy release maintenance. The energy 
loss is compensated by the external driving, that is random 
addition of the energy to the system. We do not assume that 
the driving is weak and therefore there is no complete time 
separation, so that the system does not have to be in a self- 
organized critical state, but is a driven avalanching system. 
Here we add the diffusion which proceeds independently of 
the threshold. The dynamical equations for this continuum 
burning model with diffusion can be written as follows: 



dT 
~dt 



Q(a-l) + div J+r/, 



J = giad(-Q+DT), (1) 



where time and length are properly normalized, D is the con- 
stant diffusion coefficient, and -q is random external driving. 
The energy release Q nonlinearly depends on the temperature 
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and is also history dependent (hysteresis behavior): 

Qi = kTi[9(Ti - T c ) + 9(T C - T % )9{Ti - T )#(-T,)], (2) 

It has been shown that the dynamics of the burn- 

ing model differs substantially from the sandpile dynamics. 
Active time duration is not power-law distributed which is in- 
terpreted as absence of self-organized criticality due to the 
suppression of largest avalanches. As we mentioned previ- 
ously the mean energy losses are due to direct radiation and 
not only from the boundaries. The instability criterion is local 
and does not depend on the behavior of the members. One 
of the main differences from sandpiles is that the avalanching 
activity occurs homogeneously throughout the system, while 
in a sandpile it is always enhances toward the lower end and is 
nearly zero (no active sites) near the top. The last is the reason 
for the spatially separated diffusion dominated and avalanche 
dominated regions observed in Ref. 12 and for the edge trig- 
gering of large events. 

The nonlinear equations for the burning model dynamics 
cannot be solved analytically. In order to perform numerical 
simulations we reformulate the model in the discrete form of 
one dimensional cellular automaton of L sites. At each time 
step amount of energy q enters the system at a randomly cho- 
sen position (site number i, i = 1, L) with the probability 
p. Thus the average energy input into a single site is qp and 
the total time average driving into the whole system is qpL. 
The energy release from the site i at the time step t due to the 
burning is 

Q i {t) = kT{t)[9{T{t)-T c ) 

+ 6{T C - T{t))8(T{t) - Ti)0(Qi(t - 1)] , 

and the total energy flux into site i is 

T t {t + 1) - Ti{t) = -Qi + |(1 - a)(Qi_i + Q l+1 ) 

+ D[{T i+1 - Ti) + (Ti_i - T t )\ + rj, 

(4) 

where the term with D is the regular diffusion. The role of 
the diffusion is to smooth out inhomogeneities of the tem- 
perature, thus reducing the probability of the avalanche trig- 
gering in a single site. Simplifying the explanation, let us 
assume that all sites are initially at some average tempera- 
ture. In the absence of diffusion, addition of energy to site 
randomly makes the probability of triggering an avalanche a 
rather inhomogeneous function of position. That is, if a site 
starts burning, the chances to trigger burning in its neighbor 
are typically not high, so that small size avalanches dominate. 
As the diffusion increases, diffusive redistribution of energy 
and smoothing out the temperature differences becomes pro- 
gressively faster. When diffusion is sufficiently high, the ran- 
domly added energy is redistributed among the neighbors as 
a result of this diffusion, faster than energy enters the system, 
so that temperature of a large number of sites increases simul- 
taneously (but more gradually). Once one of these sites starts 
burning the neighbors become unstable too, so that a large 



avalanche is triggered. In this way, small size avalanches be- 
come suppressed. Thus, the transition from weak diffusion 
to strong one will change the behavior of the system toward 
larges avalanches. The diffusion at which the transition starts 
may be estimated as follows. Let the average temperature of 
the sites be T av . The typical time that takes to a site to be- 
come active is r ac ~ (T c — T av )/qp. The typical diffusive 
time for the propagation of the inhomogeneity to distance L is 
Tdif ~ L 2 /D. Effects of diffusion should become noticeable 
when the time that takes to one site to increase his temperature 
and start burning becomes comparable or less than the time for 
the diffusive equilibration of the temperature with its neigh- 
bor: r ac >r dif ^D> L 2 qp/(T C - T av ), where L ~ 1. At 
this critical diffusion the smallest, single site, avalanches start 
to disappear. With the increase of diffusion the size of the 
smallest avalanche starts to increase too. Eventually, at suffi- 
ciently strong diffusion, a typical avalanche should cover most 
part of the system, causing efficient draining of the energy, so 
that the temperature drops down to the lower threshold value 
Ti . The system should relax to the initial state (the one before 
the avalanche) during the time r re i ~ (T av — T{)/qp. The 
system becomes dominated by large bursts. Burning starting 
at some site proceeds in both directions adding two sites at 
each step since all (see below) sites are near the threshold. 
A site keeps burning during the time T m i n /("l — a) since the 
flux to neighbors is balanced by the back flux from them, and 
the only energy loss is due to radiation. Thus, we expect that 
the avalanche in the (x, t) representation would have approxi- 
mately constant with in ^-direction. Each large burst drains the 
system leaving it at the temperature Tj. Next burst will start 
when the whole system is brought to the threshold T c , except 
a part near the edge of the length l e , where the energy accu- 
mulated during the relaxation time, w l e qpr re i is removed by 
the constant diffusive flux from the edge, w DTir re i, so that 
l e » DTi /qp. Respectively, we expect that the largest num- 
ber of active size during the avalanche would be » L — 2l e , 
where L is the size of the system. These estimates have been 
checked with numerical simulation on a medium size system 
of the length L = 100. The system parameters were cho- 
sen as follows: the upper critical temperature T c = 50, the 
lower critical temperature TJ = 15, the burning rate coeffi- 
cient k = 0.3, the radiative loss fraction 1 — a = 0.1, the 
amount heat added at each step q = 2 with the probability 
p = 0.005. The probability is chosen relatively high to ensure 
that the system is sufficiently active. 

We begin the comparison of the weak and strong diffusion 
with the time dependence of the avalanche size (total number 
of burning sites at each step). Figure^shows the time series 
for no diffusion D = (left) and strong diffusion D = 0.005 
(right). In the case with diffusion avalanches are less frequent 
than in the case without diffusion but each avalanche includes 
more active sites. As expected, the diffusion process causes 
smoothing of the temperature gradient so most of the sites 
have temperature around the same value and when burning 
begins in one site, it eventually sweeps almost all the system. 
Moreover, since the time of burning of almost whole system 
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FIG. 1 : Avalanche size for the case without diffusion (top) and with 
strong diffusion (bottom). 

is the same, and the relaxation time depends only on the driv- 
ing, the avalanching activity becomes quasi-periodic with the 
period given by r re ;. 

Figure [2] shows an avalanche in space and time for the 



100r 




7250 7300 7350 7400 7450 7500 

fimae 

FIG. 2: Avalanche in space and time. Width in t-direction and the 
largest number of active sites are consistent with the analytical esti- 
mates. 

strong diffusion case. The above estimates give the width of 
40 time steps and the length of inactive edge regions of about 
15, which agrees well with the simulation. 



One of the most accepted and reliable way to learn about 
the avalanching system dynamics is to study the active phase 
duration distribution. Figure[3]shows the comparison of such 




active time duration 

FIG. 3: Distribution of the active phase duration. Stars - no diffu- 
sion, diamonds - weak diffusion, and circles - strong. 

distribution for the system without diffusion (stars), with weak 
diffusion D = 0.00005 (diamonds), and with strong diffusion 
D = 0.005 (circles). It is clearly seen that weak diffusion 
starts to "eat up" the smallest durations. In the strong diffusion 
case there are already no short periods of activity but a clear 
cut maximum around t 100. 

Yet another convenient variable for comparison is the av- 
erage temperature of the system. We expect that the average 
temperature will fluctuate around some mean value in the case 
of weak diffusion. In the case of strong diffusion long periods 
of gradual increase should be followed by abrupt discharges. 
The amplitude of variations should be much larger, so that 
the maximum temperature should be quite close to the upper 
critical temperature, while the minimum temperature should 
be close to the lower critical temperature. Indeed, these ex- 
pectations are confirmed by numerical simulations, as is seen 
in Figure |3 In this figure we present the average tempera- 
ture with the Morlet wavelet transform|15ll. In contrast with 
the Fourier transform, wavelet transform allows separation of 
contributions from the narrow temperature peaks themselves 
and from the approximate periodicity of the appearance of 
these peaks themselves, provided (as expected) that the width 
of the peaks is substantially smaller than the distance between 
the peaks. For our present purposes the wavelet transform pre- 
sentation is enhanced to emphasize larger scales. As we can 
see in the case with strong diffusion, there is clear periodicity 
with the period of about 150 x 16 time steps. In the case with 
no diffusion we do not see such kind of behavior. 

This behavior we see also in the case of the number of ac- 
tive sites (Figure^ 

In conclusion, in this letter we have shown that the role 
of the diffusion process is to smooth out the temperature dis- 
tribution in the system. As a result it enhances avalanch- 
ing activity by causing large bursts at the expense of smaller 
ones. Strong diffusion even changes the mode of behavior 
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providing conditions for avalanche spreading, in the case of a 
sandpile the avalanching activity is related to the spatial sep- 
aration of the diffusion dominated and avalanche dominated 
regions. Large quasi-periodic avalanches start at the lower 
edge of the system and propagate up the slope to the inter- 
face between the two regions. In our case there is no spatial 
separation of different regions, and the spatial distribution of 
activity is not triggered at the edges. The switch in the behav- 
ior of the system is not threshold like but is gradual due to the 
suppression of the small size avalanches for strong diffusion. 
Periodicity of the bursts in the strong diffusion case is related 
to the draining of the whole system during the large avalanche 
with subsequent gradual filling, and the period is determined 
by the external driving and separation of the upper and lower 
hysteresis thresholds. 

The author is grateful to M. Gedalin for useful discussions 
and help in preparation of the manuscript. 
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FIG. 4: Average temperature as a function of time with the corre- 
sponding wavelet (Morlet) transform, without diffusion D = (two 
top panels) and with strong diffusion D — 0.005 (two bottom pan- 
els). Periodicity is clearly seen for the strong diffusion case. 



of the system, where small size avalanches completely disap- 
pear and the system enters the regime of quasi-periodic bursts 
which span most a large part of the whole system. The max- 
imum number of active sites is determined by the number of 
sites effectively drained by diffusion at the edges. The typical 
avalanche size is determined by the propagation of burning 
from the ignition point towards the edges of the active region 
and the duration of the burning of a single site with radia- 
tion losses only. The time between the two successive bursts 
depends only on the difference between the upper and lower 
critical temperatures and driving. It is important to emphasize 
the difference of the found quasi-periodic behavior from that 
described in Ref. 12. Although in both cases the basic role 
of the diffusion is to smooth out inhomogeneities and thus 
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